Investigation of g-dependent dynamical heterogeneity in a colloidal gel by x-ray 

photon correlation spectroscopy 



O 

o 

(N 

> 
O 



o 



o 
o 



> ■ 
O 

i> : 
o . 

a ■ 



I 

O 

o 



X 



V. Trapped, E. Pitard^, L. Ramos^, A. Robert^^"^ jj. Bissigi, L. Cipelletti^^* 
^ Departement de Physique, Universite de Fribourg, Chemin du Musee 3, 1700 Fribourg, Suisse 
^Laboratoire des Collotdes, Verres et Nanomateriaux (UMR CNRS-UM2 5587), 
CC26, Universite Montpellier 2, 34095 Montpellier Cedex 5, France 
'^European Synchrotron Radiation Facility - 6 rue Jules Horowitz BP 220, F-38043 Grenoble Cedex 9, France 
Present address: Stanford Linear Accelerator Center 2575 Sand Hill Road, Menlo Park, CA 94025, USA Q 

(Dated: February 6, 2008) 

We use time-resolved X-Photon Correlation Spectroscopy to investigate the slow dynamics of col- 
loidal gels made of moderately attractive carbon black particles. We show that the slow dynamics is 
temporally heterogeneous and quantify its fluctuations by measuring the variance x oi the instanta- 
neous intensity correlation function. The amplitude of dynamical fluctuations has a non-monotonic 
dependence on scattering vector q, in stark contrast with recent experiments on strongly attractive 
colloidal gels [Duri and Cipelletti, Europhys. Lett. 76, 972 (2006)]. We propose a simple scaling 
argument for the g-dependence of fluctuations in glassy systems that rationalizes these findings. 

PACS numbers: 64.70.Pf, 82.70.Gg, 82.70.Dd 



I. INTRODUCTION 

Understanding the dramatic slowing down of the dy- 
namics in systems undergoing a glass transition is one 
of the key problems in condensed matter and statisti- 
cal physics In recent years, research efforts have fo- 
cused on the role of dynamical heterogeneity: as the glass 
transition is approached, the dynamics becomes increas- 
ingly correlated in space, since rearrangements are pos- 
sible only through the cooperative motion of "clusters" 
of particles [2, 3, JJ. This cooperativity leads to strong 
temporal fluctuations of the dynamics. Indeed, because 
of dynamical correlations, the number of statistically in- 
dependent objects in the system becomes smaller than 
the number of particles, leading to enhanced fluctuations. 

Numerical simulations have tested these features on a 
wide variety of systems Experimental work, by con- 
trast, is much more scarce, because probing the dynam- 
ics with the spatial and temporal resolution needed to 
highlight their heterogeneous nature is an arduous task, 
especially for molecular glass formers 0, Q. The slow 
dynamics of colloidal systems, foams and granular mate- 
rials [5] share intriguing similarities with those of glass 
formers, including dynamical heterogeneities. These are 
experimentally more accessible than in molecular sys- 
tems, since the relevant time and length scales are larger. 
Various techniques have been used to characterize them, 
from direct space measurements d, 0, [1] to novel scat- 
tering methods that probe the temporal fluctuations of 
the intensity correlation function [l^, [HI, [13 • 

Valuable information on the physical origin of the aver- 
age dynamics is generally obtained by studying its length 
scale dependence, e.g. the dependence of the intensity 
correlation function on the magnitude of the scattering 
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vector q in scattering experiments. Similarly, one ex- 
pects to gain a better understanding of dynamical het- 
erogeneities by analyzing their behavior at different q's. 
Unfortunately, experimental and numerical or theoretical 
determinations of dynamical fluctuations as a function of 
q are still very scarce [tI, |12|, llM Q El, K [13 > leaving 
this issue an open question. 

In this paper, we investigate dynamical fluctuations 
in colloidal gels made of moderately attractive carbon 
black (CB) particles, to which a dispersant is added to 
control the strength of the (attractive) interparticle in- 
teractions. We apply, to our knowledge for the first time, 
time-resolved scattering methods to X-Photon Correla- 
tion Spectroscopy (XPCS), thereby demonstrating that 
the dynamics of the CB gels are temporally heteroge- 
neous. Dynamical fluctuations are quantified by means of 
a g-dependent dynamical susceptibility, x, similar to the 
dynamical susceptibility X4 studied in simulations [l^. 
Surprisingly, x is found to initially increase with q, but 
eventually to decrease at large scattering vectors. This 
non-monotonic behavior is in contrast with recent low- 
q measurements on diluted, strongly attractive colloidal 
gels [T^ , where dynamical fiuctuations increased linearly 
with q over one decade in scattering vector. We pro- 
pose a simple scaling argument for the g-dependence of 
dynamical fluctuations in glassy systems, which recon- 
ciles these contrasting findings and rationalizes previ- 
ously published data for granular media and glass for- 
mers [3,[il,[ii,[il,[3. 

The paper is organized as follows: in Sec.|TT]we present 
and characterize our experimental system and introduce 
shortly the time-resolved XPCS technique. Section IIIII 
reports both the average dynamics of the CB gels and 
its temporally fluctuations. Our results are discussed in 
Sec. IIVI where a simple scaling argument for the length 
scale dependence of the dynamical susceptibility is intro- 
duced. 
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II. MATERIALS AND METHODS 
A. Sample preparation and characterization 

Particle size and morphology. Our gels are made of 
CB particles of average radius R — 180 nm, suspended 
in mineral oil at an effective volume fraction ip f« 6%, as 
determined by viscosity measurements (see below). The 
particles have a conveniently high scattering cross sec- 
tion for X-rays. They are fractal aggregates made of 
permanently fused primary particles (fractal dimension 
df = 2.2±0.1) [l^. The diameter of the primary parti- 
cles ranges from 20 to 40 nm, as determined by electron 
microscopy. However, the smallest units in our samples 
are effectively the CB particles themselves and not the 
primary particles, since the CB particles can not be bro- 
ken, neither by thermal fluctuations, nor by adding a 
dispersant or by applying a large shear, e.g. during rhe- 
ology tests or sonication. The size and polydispersity of 
the CB particles were determined by applying the cumu- 
lant analysis described in Rcf. 19] to intensity correlation 
functions measured by dynamic light scattering (DLS) at 
q = 6.01 /im~^ (scattering angle 6* = 20 deg). A very 
diluted ((^ « 2 X 10~^) and fully dispersed sample of CB- 
particles in mineral oil was prepared for the DLS mea- 
surements. For spherical, monodisperse particles, one 
expects (72 ~ 1 to relax exponentially; deviations due to 
shape and/or size polydispersity may be quantified by the 
ratio K2/K1, where Ki is the i—th coefficient of a cumulant 
expansion: ln{[g2(''") — 1]°'^} — kq — kit+K2t'^ /2 + ... p^ . 

Figure [T] shows ln{[(72(T) — l]'^'^} vs t for our DLS 
measurements: the data are very close to a straight line, 
indicating a nearly exponential decay and thus suggest- 
ing that the sample polydispersity must be moderate. By 
generalizing the arguments of Ref. [l^ to particles with 
a fractal morphology, one finds that in the low-g limit 
{qR < 1) the first two cumulants are related to the aver- 
age translational diffusion coefficient, D, and its relative 

variance, a'^ 
sions: 



{D'^/D — 1), by the following expres- 



Ki = Dq^ 



(1) 
(2) 



with D = kBT/{6TTf]R,pp), i?app = i?2'if/i?2rff-i, T the 
temperature and fee the Boltzman's constant. Here, 
i?" — J dRP{R)R'"' is the n-th moment of the normalized 
number distribution of particle radii, P{R). For the data 
shown in Fig. 1, we find aj-, — 2k2/h\ = 0.047, confirm- 
ing that the polydispersity is moderate. The average par- 
ticle radius and its relative variance, cr|j — {R'^/R — 1), 
may be obtained from and i?app provided that P{R) 
is known. Since for moderate polydispersity the exact 
shape of P{R) has little influence on the final result, we 
choose a generalized exponential (or Schulz) distribution, 
for which calculations can be performed analytically . 
By taking for simplicity d{ = 2 (very close to the value 




FIG. 1: Second-order cumulant fit (line) of the intensity cor- 
relation function (open circles) measured by dynamic light 
scattering at q = 6.01 fim"^ for a diluted suspension of CB 
particles. 



2.2 of our particles), one finds 



i?app = RVR^ = i?(l + 34) 



i?2 fl4 



1 + 34 

1 + 24 



(3) 
(4) 



yielding for our CB particles R = 180 nm and = 0.23. 

Determination of the volume fraction The particle 
concentration of the two CB gels studied here is 2% 
w/w. The effective volume fraction corresponding to 
this weight fraction was determined by measuring the 
viscosity of a suspension where the CB particles were 
fully dispersed (1.6% w/w dispersant). The effective hy- 
drodynamic volume fraction, (p, is estimated using [20| 
Voo/tJ' = (1^',^'/0.71)^^, where /i is the shear-independent 
viscosity of the mineral oil and rj^ is the high shear rate 
viscosity of the suspension. The value (p ~ 6% thus ob- 
tained was found to be consistent with that obtained by 
measuring the low-shear viscosity of the suspension, ryo, 
and using [2d| t/o/a* = (1 - ^fi/O.eS)-^. 

Estimate of the depth of the particle-particle interac- 
tion potential. The attraction between CB particles is 
controlled by the amount of added dispersant. The depth 
of the well of the particle-particle interaction potential, 
[/, may be estimated by studying the mechanical re- 
sponse of the CB gels. As discussed in ref. ^l|, the relax- 
ation spectrum (as measured by oscillatory rheology) of 
CB suspensions and its dependence on volume fraction is 
essentially identical to that obtained for model systems of 
nearly monodisperse spherical particles that interact via 
a short-ranged potential induced by the depletion mech- 
anism. Therefore, we use the non-equilibrium state dia- 
gram of the depletion systems with known depth of the 
interaction potential to evaluate U for our CB gels. The 
diagram is shown in Fig. 3 of Ref. [2l| : the control pa- 
rameters are ip and U; a line U = Uc{f) separates the 
systems with solid-like mechanical behavior (upper part 
of the diagram) to those with a liquid-like behavior (lower 
part). For the CB gel with 0.2% dispersant, we find that 
the elastic modulus is very low {Gp = 9 x 10""* Pa), 
which indicates that this system corresponds to a state 
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close to the fluid-solid boundary. Figure 3 of Ref. [2l| 
indicates that for short ranged potentials the boundary 
corresponds to U ^ Uc ~ 12kBT at ip — 6%. Thus, we 
assign U ~ 12fcBT to the CB gel at (p = 6% and with 
0.2% dispersant. 

For the CB-system with 0% dispersant, we take advan- 
tage of the dependence of the elastic modulus on disper- 
sant concentration, which was found to exhibit a critical- 
like behavior, Gp = Go{U/Uc — 1)", with h' = 3.9 and 
Go « 1 Pa We measure Gp = 4.5 Pa for the CB gel 
without dispersant, yielding U w 2.5Uc « 30 fceT. 

Sample preparation. The CB suspensions are prepared 
from a stock suspension of 4% w/w carbon black (Cabot 
Vulcan XC72) in light mineral oil (Aldrich). The stock 
solution is made by first mixing the carbon black powder 
into the mineral oil using a standard household mixer. 
Subsequently, it is thoroughly sonicated for 1 hour using 
an ultrasound device (Hielscher UP200H) operating at a 
power of 120 W and a frequency of 24 kHz. To avoid 
overheating, the suspension is placed in an ice bath and 
sonication is run in cycles of 1 s. The stock solution is 
then left to equilibrate for 2 days, before it is used to 
prepare the final samples. Before diluting the sample 
to a final concentration of 2% w/w CB (corresponding 
to (fi ^ 6%), the stock suspension is again thoroughly 
mixed; the final samples are then prepared by adding 
either pure mineral oil or a solution of dispersant in min- 
eral oil; they are mixed and left to equilibrate for 2 days. 
Prior to the XPCS experiments, the samples are again 
thoroughly mixed and subsequently injected in a 1 mm 
(inner diameter) cylindrical capillary. The capillary is 
sealed and placed in a water bath in which we introduce 
our ultrasound device, thereby submitting the sample to 
a final indirect sonication step (120 W and 24 kHz) for 
about half an hour. We define the age of sample, t^, 
as being the time elapsed since this final sonication step. 
Note that sonication breaks up any aggregates that may 
have been formed between the CB particles, without af- 
fecting the integrity of the CB particles themselves. 



B. Time-resolved XPCS measurements 

Average dynamics. The dynamics of the CB gels are 
investigated by means of XPCS performed at the IDIO 
Troika beamline at the European Synchrotron Radiation 
Facility (ESRF), with X-rays of wavelength A = 1.55 
A. The scattering volume has a cylindrical shape, with 
diameter 12 fim and length 1 mm (along the horizon- 
tal direction), defined by the beam size and the capil- 
lary diameter, respectively. The scattered intensity is 
recorded by a charge-coupled device (CCD) camera, cov- 
ering about one decade in scattering vector: 15 ^m~^ < 
q < 150 /im~^, corresponding to distances comparable 
to or smaller than the CB particle size. XPCS measure- 
ments start at age t^, = 2 min: we find that the dynamics 
initially slows down, as observed for many glassy sys- 
tems [5|] , but that after a few hours a stationary state is 



attained. All measurements presented in this work start 
at tjju = 5 hours, well in the stationary regime, and last 
typically up to 3 hours. Mechanical and beam instabil- 
ities are a concern when measuring very slow dynamics, 
especially in XPCS experiments. By using a static scat- 
terer (Vycor glass), we have checked that the setup was 
stable up to T « 3000 sec, longer than the relevant time 
scales in our experiments. 

In order to measure the average dynamics and 
its fluctuations, we use the time-resolved correlation 
(TRC) method [9|. The instantaneous degree of 
correlation between X photons scattered at time t 
and t + r is measured according to ci{q,t,T) — 

{Ipit)Ip{t + T))^ I ((/p(0), (/p(i + r)) J - 1, where Ip(t) 

is the scattered intensity at pixel p and time t and (• • •)^ 
is an average over a ring of pixels corresponding to ap- 
proximately the same magnitude of q but different az- 
imuthal orientations. The intensity correlation function 
is (72(9, ■'') ~ 1 = cj, where indicates a time average. 
gi is related to the dynamic structure factor /(g, r) by 
gi — t — /3/^, where /3 < 1 is a positive instrumental 
constant. 

Fluctuations of the dynamics. To quantify dynamical 
heterogeneity we calculate xi''': q)'"'^^^'' = var(c/(g, r)), 
the temporal variance of the instantaneous degree of cor- 
relation, ci{q,t,r) [l^l, which we correct by subtract- 
ing the contribution of measurement noise. The correc- 
tion procedure is described in detail in Ref. [2^; here 
we simply recall its main features. The experimentally 
measured degree of correlation is affected by a statistical 
noise stemming from the finite number of pixels, rip, over 
which this quantity is averaged. Therefore, one has [2^ 

x(t, 9) = X^^'^pHt, g)-var[n(g, t, r)] = x^°'^P^(r, <z)-G/np , 

with x^'^^^\'''^q) = var[c/(g, i, r)], x(t, q) the desired 
noise-free dynamical susceptibility, n{q, t, r) the statis- 
tical noise, and C a positive constant. The last equality 
stems from the central limit theorem. In order to obtain 
x(t, g), for each q we calculate x^°^'^''{'^^ q) by processing 
the same set of images with different choices of rip (by 
using only 1 pixel every 1, 2, 4, 8,... available pixels). By 
extrapolating the X^^^^Ht^q) I/^^p data to l/rip — 0, 
we obtain the desired noise- free x{t, q)- To estimate the 
uncertainty associated with this procedure, we inspect 
the data for r much larger than the slowest relaxation of 
.92 — 1, where one expects x(t, g) = [l3, UHl- In this 
limit, we indeed find x = (see Figs. 2b and 3b below) 
within an experimental uncertainty, a^{q), which is es- 
timated by calculating the standard deviation of x(''"j'Z) 
as a function of r, at large r. The values of cFxil) thus 
obtained are taken as an estimate of the error bars on x 
at all delays, and will be used in particularly in Fig. 5. 
The average relative uncertainty is 6% for the CB gel 
with U ^ 12 k-oT and less than 1% for the gel with 
t/ « 30 fceT. 
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FIG. 2: (Color online) Top panel, main figure: intensity corre- 
lation functions at various q's for the CB gel with (7 ~ 12 fceT 
(data are corrected for stray light). From top to bottom, q 
increases from 17.1 jj.m~^ to 92.0 jj.m~^ (see labels in bot- 
tom panel). Inset: g-dependence of ri (open squares) and T2 
(solid squares). The line is a power law fit to ri(g) yielding 
an exponent —0.91 ± 0.1. Bottom panel: r dependence of the 
dynamical susceptibility x for the q vectors corresponding to 
those shown in the top panel. For the sake of clarity, not all 
available data have been plotted in the two main panels. 



III. RESULTS 

Figure 2a shows 52 — 1 for the CB gel with U ~ 12 k^T, 
at various q. The intensity correlation function exhibits 
a two-step decay, whose physical origin we shall dis- 
cuss later. Figure [^b displays the r-dependent ampli- 
tude of the noise-corrected fluctuations of the dynamics, 
x(t, g); at all g's, x exhibits a peaked shape, strongly 
reminiscent of that observed in a variety of glassy sys- 
tems 0, S [H [H [H El [H, m, Hi. The peak of the 
dynamical susceptibility, %*, occurs at a time delay T*{q) 
of the same order of magnitude of the decay time associ- 
ated with the initial relaxation of g2 — 1. Additionally, a 
shoulder in x(t) is observed on the time scale of the final 
relaxation of 32 — 1- 

To gain insight on the physical origin of the slow dy- 
namics, it would be desirable to quantify the decay of 52 
by fitting its relaxation by an appropriate function, such 
as the combination of two (stretched) exponential func- 
tions. However, due to the limited image acquisition rate, 
the initial decay of g2 — 1 is only partially captured in 
our experiments, especially at the largest q vectors. As a 
consequence, a direct fit of the initial relaxation is not re- 
liable. As an alternative procedure, we use a scaling anal- 



ysis similar to that adopted in several photon correlation 
spectroscopy works (see, e.g., [H, [23, HI], as shown for 
the gel with J7 « 12 k-BT in Fig. 3a. We first scale the lag- 
time axis by t* (g) , the lag corresponding to the peak of 
the dynamical susceptibility x- This choice is motivated 
by the fact that, quite generally, for glassy systems t* is 
of the same order of magnitude of and proportional to the 
system relaxation time (see e.g. [1, 129]). Furthermore, 
our X data exhibit a remarkable scaling behavior, as 
shown in Fig. 3b. This allows us to determine r* at all q 
with an uncertainty smaller than the temporal resolution 
of the CCD camera (2 sec). We finalize our scaling proce- 
dure of the correlation function by normalizing 32 — 1 with 
a g-dependent amplitude, A{q), which leads to a collapse 
of all our (/-dependent data in the lag-time-range of the 
initial decay. To determine the shape of the initial decay, 
we assume that the data can be modeled by a stretched 
exponential decay, ^ exp[— (T/ri)Pi], and search for the 
stretching exponent pi that best linearizes our data plot- 
ted as ln[(ji2('?, t) — 1] — InA(g) versus (t/t*)^\ where we 
find an optimum value of pi — 1.2 ± 0.2. As shown in 
Fig. 3a, a very good scaling onto a straight line is ob- 
tained for T < r*; at larger lags the data curve upwards, 
because the initial decay of 32 is followed by a (tilted) 
plateau, as seen in Fig. 2a. The same scaling proce- 
dure is also used to characterize the initial decay for the 
gel with f7 « 30 k^T, for which we find pi = 1.1 ± 0.2 
(data not shown). Figure 4 shows some representative 
intensity correlation functions for both gels (symbols), 
together with the fits issued from the scaling analysis 
(lines). A full relaxation of 52 is observed only for the 
gel with U K, 12 k-gT, whose final decay is well approxi- 
mated by a single exponential as shown in Fig. 4b. The 
fact that for the gel with [/ « 30 k-oT no final relaxation 
of §2 is observed in the accessible time scale is most prob- 
ably due to the deeper interparticle potential well that 
makes particle displacements more difficult. 

The g-dependence of ti , the characteristic time of the 
initial relaxation of 52 ^ 1, is shown in the inset of Fig. [2^ 
for the gel with U ^ 12 k-gT. The data can be modeled 
by a power law ri ~ g-o.9i±o.i^ shown in the inset of 
fig. [2^ as a continuous line (similar results are also ob- 
tained for the gel with U ~ 30 fceT). We can exclude 
that the initial relaxation of §2 is due to thermally in- 
duced fluctuations of the gel branches [l^l, which we ex- 
pect to exhibit characteristic times at least three orders 
of magnitude smaller than ti. Indeed, both the nearly 
g^^-dependence of ri and pi > 1 suggest that the dy- 
namics is determined by stress-induced rearrangements 
similar to those of other soft glassy materials (see 5, 2^ 
and references therein). The final decay of the correla- 
tion function is approximately exponential, with a decay 
time T2 that is (/-independent. This behavior most prob- 
ably stems from random rare rearrangements occurring 
when bonds are broken, leading to particle displacements 
larger than 1/q. Indeed, in this case we expect one sin- 
gle rearrangement to be sufficient to fully decorrelate the 
contribution of the displaced scatterers to g2 — 1- The 
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FIG. 3: (Color online) Scaling analysis of the average dynam- 
ics and of its fluctuations, for the CB gel with [/ ~ 12 k^T 
(same data as in Fig. 2; here data at all the available q are 
plotted, except for the smallest and the largest, for which the 
data are more noisy), a) Scaling plot of the initial decay of 
g2 — \. b) Scaling-plot of the dynamical susceptibility using 
reduced variables x/x* s^nd t/t*. x* s^nd r* are the height 
and the position of the peak of the dynamical susceptibility, 
respectively, (same symbols as in a)). 



only time scale for the final relaxation of g2 is then the 
average time between rearrangements, regardless of q. 
Assuming uncorrelated events with Poissonian statistics, 
one expects an exponential decay of 52, in agreement with 
our measurements. As mentioned before, for the CB gel 
with no dispersant {U « 30 k-oT) no final relaxation is 
observed within our experimentally accessible time win- 
dow; this is most likely due to a decrease in the rate of 
bond breaking as U increases. 

We investigate the length scale dependence of dynam- 
ical heterogeneity by measuring x*{q)i the height of the 
peak of the dynamical susceptibility, shown in Fig. 5. 
For both gels, we find that x*{l) first increases with g, 
reaches a maximum value a,t q = q* and then decreases 
at larger scattering vectors. Moreover, we find that the 
dynamical fluctuations are more pronounced for the gel 
with the deeper interparticle potential well. This can be 
intuitively understood as the size of the regions that re- 
arrange cooperatively is presumably larger in gels with 
stronger interparticle interactions, leading to larger dy- 
namical fluctuations. Additionally, it is conceivable that 
the relaxation process itself is more heterogeneous in time 
for stronger gels, where internal stress build-up and re- 
lease is likely to be more important, further contribut- 




FIG. 4: (Color online) a): Intensity correlation functions rep- 
resentative of the behavior at small, intermediate and large q 
vectors for the CB gel with 0.2% dispersant (symbols; from 
top to bottom, q — 17.1,39.6, and 80.8 /im~^). The lines 
are fits to the initial decay of — 1 issued from the scaling 
analysis described in the text, b): zoom on the final decay of 
(72 — 1 for the same data as in a), together with the large-r fits 
(lines) . c) : g2 — \vsT for the CB gel with no dispersant (sym- 
bols; from top to bottom, q = 26.5,61.2, and 141.5 /im"^). 
In a) and c) the fits are represented as solid lines for r < ri 
and continued as dashed lines beyond the fitting interval. 
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FIG. 5: Peak of the dynamic susceptibility, x*, 9 for a CB 
gel with 0.2% dispersant (squares) and no dispersant (circles). 
Error bars are calculated as explained in sec. [Ill Inset: same 
data normalized by the squared amplitude of (72 — 1. 



ing to enhanced fluctuations. As outlined in ref. [1^ . x 
should in principle be normalized by the (g-dependent) 
squared amplitude of the relaxation of 172 — 1, to allow 
for an explicit comparison of data obtained at different 
g's. [11]. Here, such a correction is affected by a large 
uncertainty, since the short-time behavior of (72 is barely 
accessible to the CCD. The height of the peak of the 
variance normalized using the amplitude estimated via 
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the scaling procedure, x*/A^, is shown in the inset of 
Fig. [5l Although the data are more noisy, the peaked 
shape observed in the main plot is preserved. 



IV. DISCUSSION 

The q dependence of x* found here is very different 
from that measured in low-g light scattering experiments 
on strongly attractive gels made of spherical polystyrene 
particles [12l |. where x* ~ 1- One may wonder whether 
this discrepancy stems from a difference in the particle 
morphology or polydispersity. However, CB suspensions 
have been found to exhibit the same relaxation spectrum 
(as measured by rheology) as model systems consisting 
of spherical, nearly monodisperse colloids interacting via 
short-ranged depletion forces [2l|. Thus, the dynami- 
cal properties reported here are likely to be generic for 
colloidal gels with moderate attractive interactions. In- 
stead, we propose that the different g-dependence of x* 
arises from the different q and U ranges probed in this 
work compared to p^ . In this section, we introduce a 
simple — yet general — scaling argument for the length- 
scale dependence of dynamical fluctuations in glassy sys- 
tems that reconciles these contrasting observations. We 
assume the dynamics to be due to random rearrangement 
events each affecting a "blob" of volume Vh- The fluc- 
tuations of c/ then stem from fluctuations of iVtot, the 
total number of events needed to decorrelate the scat- 
tered light. In determining TVtot, one has to take into 
account that one single event may not displace particles 
far enough to fully suppress the local contribution to the 
scattered light; additionally, one single event may affect 
only a portion of the whole scattering volume. Thus, 
^^tot ~ ^biob^Gv, where A'^biob ~ Kic/Vb is the number of 
dynamically correlated blobs of volume Vb contained in 
the scattering volume Vsc , and iVov is the number of rear- 
rangement events that are needed, at any given location, 
to relax the local contribution to the correlation function, 
i.e. the number of events on the time scale Ti- of the sys- 
tem's relaxation. The inset of Fig. [S] is a schematic rep- 
resentation of this concept for iVbiob ~ 50 and A^ov = 2: 
in this case, 92 — 1 decays to zero when at least two rear- 
rangement events (symbolized by a black and a gray cir- 
cle) have occurred at every location in the system. Given 
the random nature of the events, A^tot fluctuates, and so 
does c/, with var(c/) var(A'tot)- According to the cen- 
tral limit theorem, the relative variance of A^tot — and 
thus X — is expected to scale as N^l ^ (A^biobA^cv)"^- 

We stress that in general both A'biob and A'ov are q- 
dependent quantities. Indeed, A^cv decreases as q in- 
creases, because fewer events are needed to displace 
the particles over the smaller distances corresponding 
to larger scattering vectors. More specifically, in our 
model the average number of events is proportional to 
time and thus we expect A'ov ~ Tj. ~ 9 For uncor- 
related particle displacements due to successive events 
( "Brownian-like" rearrangements), a = 2. By contrast. 




log^ 

FIG. 6: (Color online) Inset: schematic representation of the 
rearrangement events within the scattering volume. Succes- 
sive events at the same location are indicated by circles of 
different grey level. Main figure: qualitative double logarith- 
mic plot of the proposed g-dependence of dynamical fluctua- 
tions. Left axis: x* (continuous line, blue); right axis: A*'^"^^ 
(dashed line, red), and A^,~j^i, ~ Vb (dotted line, black). The 
slope Q varies between 1 and 2, depending on the nature of 
the dynamics, as discussed in the text. 



a = 1 if the displacement direction persists over several 
events ("ballistic-like" rearrangements), as found for our 
CB gels and other systems with internal stress-driven dy- 
namics d, [13, nil . Note that at very large q N^v should 
saturate to one, when the particles' displacement due to 
one single event exceeds l/g, the length scale probed in a 
scattering experiment. By contrast, there are no a priori 
prescriptions on the g-dependence of Abiob- However, it 
is reasonable to assume that the latter be an increasing 
function of q, since in glassy systems rearrangements in- 
volving displacements over large distances (probed at low 
q) are likely to require the highly cooperative motion of 
many particles, whereas smaller displacements (probed 
at high q) may be achieved independently by clusters 
containing just a few particles. Thus, for q ^ 0, A^biob 
is expected to saturate at a lower bound, which may be 
as low as 1 when Vb is larger than the scattering vol- 
ume, as observed in 12]. In the opposite limit, 5 — *■ 00, 
A'biob should saturate to A^p, the number of particles in 
the system. 

The main panel of Fig. [5] illustrates schematically the 
g-dependence of A^biob ^cv^ (right axis) and that 
of X* ~_(A^biobA^cv)~^ (left axis). At low q, x* - 
since A^i^^,^ saturates to one whereas A^,,"^ ~ q°'. This 
is the regime observed for the strongly attractive gels of 
ref. [1|, for which A^biob = 1 and x* ~ Mq)-^ ^ 

q (a. = 1). By contrast, x* tends to a constant value at 
very large g, because N^^^ saturates to N~^, while N^^ 
saturates to one. The behavior of the dynamical suscep- 
tibility for intermediate g's depends on the detailed inter- 
play between A^biob and A"cv In Fig. [S] we have sketched 
the case where A^i^^b decreases faster than A^^"^ grows, 
yielding a peaked shape of x* ■, as observed in the exper- 
iments presented here. 

We expect these scaling arguments to hold also for 
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other glassy systems exhibiting dynamical fluctuations. 
Indeed, a growing trend in the low q regime, with a — 2, 
is predicted for dynamically facilitated models, for which 
X* saturates in the opposite limit q ^ oo flS*] . For a 2-D 
granular system [2] and in simulations of a Lcnnard-Jones 
glass former and a short-range attractive glass p^ . 
a non-monotonic behavior similar to that of the CB gels 
has been reported, with x* initially growing with q but 
then decreasing after going through a peak. The length 
scale corresponding to this peak is comparable to the 
inter-particle distance for repulsive systems, while it is 
shifted to smaller values, corresponding to the width of 
the interparticle potential well, for attractive glasses l^lSj] . 
Finally, the high-q regime has been explored for hard 
spheres within the mode coupling theory in where 
the amplitude of x^, a lower bound for Xii was reported 
to decrease with q around the peak of the static structure 
factor, in agreement with the behavior predicted by our 
scaling argument at large q. 

Simulations of glass formers and experiments on 2D 
granular materials suggest that x* is maximum on the 
length scale of the particle size or that of the interparti- 
cle bond. Colloidal gels present an additional characteris- 
tic length, the size of the (fractal) clusters that compose 
them. For the gels of ref. [12], x* was shown to grow 
with q for length scales intermediate between the cluster 
and the particle size. For the CB gels studied here, x* 
peaks at q* « 20 — 50 ^ni^^ (see Fig. [5]), corresponding 
to a length scale A = 2TT/q* « 125 — 300 nm, compara- 
ble to the particle size. Collectively, these observations 
suggest that the crossover length scale for dynamical het- 
erogeneity in colloidal gels is of the order of or smaller 
than the particle size, similarly to molecular glass form- 
ers and granular materials, rather than the cluster size. 
Interestingly, in our gels A shifts towards smaller val- 
ues {q* increases) when the particles are more sticky, the 
same trend as that reported in Ref. when going from 



a nearly hard-sphere system to an attractive one. Al- 
though the exact value of q* most likely depends on the 
detailed shape of the interparticle potential, it is intrigu- 
ing to note that the values found here {q*R « 3.6 for 
U w 12fcBT and q*R « 9 for U « 30fcBT, respectively) 
are comparable to those shown in [15] for a nearly-hard- 
sphere and an attractive system {q* R ~ 2.5 and q*R w 6, 
respectively) . 

In conclusion, we have shown that the dynamics of 
CB gels is temporally heterogeneous. Dynamical fluctu- 
ations increase with q, peak around the inverse particle 
size, and decrease at larger scattering vectors. This be- 
havior and that of other systems can be rationalized by a 
simple scaling argument, providing a general framework 
for understanding temporal fluctuations of the dynam- 
ics in glassy systems. Additionally, our measurements 
demonstrate that XPCS may be used to obtain quanti- 
tative information not only on the average dynamics, but 
also on its heterogenous behavior. This opens a new way 
to investigate dynamical heterogeneity in a wide variety 
of materials, possibly including molecular glass formers, 
whose characteristic length scales match those probed 
by XPCS, provided that a high enough signal can be col- 
lected. 
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